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Abstract. A highly nonlinear, fourth-order wave equation that models the continuum 
theory of long wavelength pulses in weakly compressed, discrete, homogeneous chains with 
a general power-law contact interaction is studied. Eor this wave equation, all solitary wave 
solutions and all nonlinear periodic wave solutions, along with all conservation laws, are 
derived. The solutions are explicitly parameterized in terms of the asymptotic value of the 
wave amplitude in the case of solitary waves and the peak of the wave amplitude in the case 
of nonlinear periodic waves. All cases in which the solution expressions can be stated in an 
explicit analytic form using elementary functions are worked out. In these cases, explicit 
expressions for the total energy and total momentum for all solutions are obtained as well. 
The derivation of the solutions uses the conservation laws combined with an energy analysis 
argument to reduce the wave equation directly to a separable first-order differential equation 
which determines the wave amplitude in terms of the travelling wave variable. This method 
can be applied more generally to other highly nonlinear wave equations. 
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1. Introduction 

Travelling waves are an interesting feature of many nonlinear phenomena and have numer¬ 
ous important physical applications. One area that has attracted much recent attention is 
the study of compression pulses in one-dimensional systems of discrete macroscopic particles 
that interact by a power-law contact potential [TH5U]. When such a discrete system is pre¬ 
compressed, then propagating solitary waves and nonlinear periodic waves can be produced 
through suitable tuning of the initial conditions [T8l[25l[30l|3Tl[35]. Pulses whose wavelength 
is much larger than the particle radius, but much smaller than the length of the chain, 
can be modeled by a continuum partial differential equation (PDE). Continuum models are 
important because they allow comparison with experimental Endings and give insights into 
general features of compression waves in the underlying discrete system. 

The simplest model describes a homogeneous chain of A > 2 spherical particles with 
mass m such that the interaction contact potential is given by the Hertz law [ID] V = 
where a is a constant which depends on the material properties of the particles, and b = 
bo + Ui — Mj+i > 0 is the (dynamical) overlap distance between adjacent particles, given in 
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terms of the relative particle positions Uj, j = 1,... ,N, and the initial overlap ^o- This model 
has an interesting generalization in which the exponent in the Hertz potential is replaced 
by an arbitrary power, V = where k > 1 for typical physical applications jHJ-HS]. 

When the initial overlap is small compared to the particle displacements, the chain is said 
to be weakly compressed. In this case, the continnnm PDE for the relative displacement 
u = —u{x,t) of the particles is given by the highly nonlinear, fonrth-order wave eqnation 

[smsusaEs 

C Uft Uxx T '^xx T ^xx'^xxx T '^'^x ^xxxx (i-T) 

with 

C = .yk{k + l)a(2i?)^+Vm, (1.2) 

a = l{k-l){k-2)R'^, /3 = l{k-l)R^, 7 = (1.3) 

where R is the radius of the particles. Note, in these expressions fll.2p - fll.3p . the constants 
a, (3, 7 have units of length-squared, and the constant c has units of speed, while it has units 
of length. 

An important remark is that the unitless strain 

V = Ux = —Ux > 0 (1.4) 

rather than the displacement it = —u, will be the variable whose dynamics will physically 
describe a propagating wave. The strain is positive because the underlying discrete dynamical 
system has no restoring forces. For solitary waves, the strain will have the asymptotic 
behaviour 

n —)■ Uo > 0 as a; —)■ ±cxd (1.5) 

where Vq is a constant background strain. The boundary conditions on the strain for non¬ 
linear periodic waves consist of 

v{x,t) = v{x + LR + T), Vx{xR) = Vx{x + LR + T) (1.6) 

where L is the wavelength and T is the period. 

Existence of solitary wave solutions for the highly nonlinear wave equation fll.lj) has been 
proved recently [3^. However, an exact analytical expression for these solutions has not 
appeared in the literature to-date. Moreover, the only known exact solution is a special 
periodic wave [32] that contains no free parameters. 

In the present paper, we derive exact analytical expressions for all solitary wave solutions 
and all nonlinear periodic solutions of this wave equation fll.lj) . We also obtain analytic 
expressions for their total energy and total momentum. The solutions are parameterized in 
terms of the asymptotic value of the wave amplitude in the case of solitary waves and the 
peak of the wave amplitude in the case of nonlinear periodic waves. 

Our method involves the following steps. First, we hnd the conservation laws admitted 
by the wave equation fll.ll) . Next, we use these conservation laws to derive a hrst-order 
ordinary differential equation (ODE) for all travelling waves, without having to do any 
explicit integrations. This reduction contains two free constants which correspond to first 
integrals related to the energy flux and longitudinal-momentum flux of a travelling wave. 
Third, we employ an energy analysis argument to classify the types of travelling waves. 
This argument shows that there are three types of solutions: (i) solitary waves with an 
exponentially decaying tail, which satisfy the asymptotic condition fll.51) : (ii) periodic waves, 
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which satisfy the periodicity condition fll.hp : (iii) cusped waves, both periodic and solitary; 
(iv) nodal waves, in which the wave amplitude vanishes at isolated points. Fourth, we solve 
the ODE in each of these different solution cases by giving an explicit integral expression that 
determines the wave amplitude in terms of the travelling wave variable. Last, we show how 
to use the ODE along with its first integrals to obtain the total energy and total momentum 
carried by these solutions, without the need to evaluate complicated integrals involving the 
explicit form of the solutions. 

The paper is organized as follows. 

In section [2l we first present the Lagrangian and Hamiltonian formulations of the highly 
nonlinear wave equation fll.ip . and we next derive its conservation laws and corresponding 
conserved quantities by a direct method which does not require the use of Noether’s theorem. 
This material also contains some general results which are applicable to the study of any 
fourth-order wave equation. 

In section [3] we apply the steps in our method to the wave equation fll.ip for an arbitrary 
exponent k > 1. This yields explicit integral expressions for solitary waves, nonlinear periodic 
waves, cusped solitary waves, cusped periodic waves, and nodal waves, which are the main 
results of the paper. We also derive conserved integral expressions for the total energy and 
total momentum carried by these different kinds of travelling waves. 

In section |U we systematically determine all cases of the power-law exponent k > 1 for 
which each of the different solution integrals can be evaluated in an explicit analytical form 
using elementary functions. This analysis yields two explicit solitary waves, for k = 2,3] two 
explicit periodic waves, one for k = 2 with zero energy-flux and non-negative longitudinal 
momentum-flux, and the other for k = 3 with non-negative energy-flux and zero longitudinal 
momentum-flux; three explicit cusped waves, two being cusped solitary waves for k = 2,3, 
and the third being a cusped periodic wave for k = 3 with zero longitudinal momentum-flux; 
two explicit nodal waves, one for k = 2 with non-negative longitudinal momentum-flux, 
and the other for k > 1 arbitrary with zero longitudinal momentum-flux. We give explicit 
expressions for the total energy and total momentum for each of these solutions. 

Finally, we make some concluding remarks in section |5l 


2. Preliminaries 

For the highly nonlinear wave equation fll.ip . we will first present its Lagrangian formu¬ 
lation which we will then use in interpreting its conservation laws. 


2.1. Variational formulation. We start by considering a general fourth-order wave equa¬ 
tion Utt = F{u,Ux,Uxx,Uxxx, Uxxxx)- Au equation of this form has a Lagrangian formulation 
if it arises from the extremals of a Lagrangian function L = -|- Lq{u, Ux, Uxx, Uxxx, Uxxxx) 

with respect to u. A Lagrangian function is not unique, because adding any total derivative 
expression to it does not change its extremals. This freedom can be used to find a Lagrangian 
of lowest possible order. 

The extremals of L are obtained by setting its variational derivative to be zero, where the 
variational derivative (Euler-Lagrange operator) is defined by [351IT6] 
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which has the property that it annihilates any total derivative expression. Here D denotes 
a total derivative (acting by the chain rule), namely 
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Hence, a Lagrangian formulation then consists of finding a function Lq for which 


f) = utt- F = = utt - (2.3) 

ou du 

It is straightforward to show that, since F is fourth-order, if Lq contains any terms of higher 
than second order then these terms must take the form of a total derivative, as otherwise 
they will produce terms of higher than fourth order in the variational derivative dL^/du. 
Consequently, there is no loss of generality if we require Lq to be only second order. 

Thus, a Lagrangian formulation fl2.3p will exist for Uu = F{u,Ux,UxxjUxxx^Uxxxx) if (and 
only if) we can find a function Lq{u,Ux,Uxx) that satisfies 


SLq dLr 


"dux "duxx 


= F. 


Su du 

The highly nonlinear wave equation fll.ip has a simple Lagrangian function 

r = Pt + ro(fir.iii). io = ■ 


(2.4) 

(2.5) 


This Lagrangian has the lowest possible order. By comparison, in Ref. [18] a more compli¬ 
cated higher-order Lagrangian was derived. It is easy to check that the difference between 
that Lagrangian and the one given here is just a total x-derivative expression. 

Next we show how to convert this Lagrangian formulation into an equivalent Hamilton¬ 
ian formulation using the strain variable (II.4p . as this will be helpful later for physically 
interpreting the conservation laws admitted by the highly nonlinear wave equation fll.ip . 

We consider a general Lagrangian wave equation of the form fl2.3p . This equation can be 
written equivalently as a first-order system in time derivatives if we introduce variables for 
the first-order derivatives of u: 


Ut = w, Ux = —V. 


( 2 . 6 ) 


Then we have Wt = Uu = F = dLo/du, and Vt = —Utx = —Wx- The variational derivative 
dLo/du on the righthand side of the Wt equation can be expressed in terms of v through the 
identity [IS] 


(2.7) 


L = dL, 

du Sv 

Using this identity, we now have Wt = Dx{5Lq/ dv). The righthand side of the Vt equation 
can be expressed in a similar form by noting Wx = Dx{d{^w‘^)/dw), so then we have Vt = 
— Dx{S{\w‘^)/dw). Thus, the original wave equation becomes 


where 
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H dx, 


H = \vp - L„ 


( 2 . 8 ) 
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is the Hamiltonian, and where 


( 2 . 10 ) 


V = 




is a skew-adjoint matrix operator. (Note the choice of sign in P corresponds to a positive sign 
for the kinetic term in H). There is a non-canonical Poisson bracket associated 

to this formulation. For any two functionals T = F dx and Q = G dx, given in 
terms of functions F{v, w, Vx, Wx, Vxx, Wxx, ■ ■ •) and G{v, w, Vx, Wx, Vxx, Wxx, ■ ■ ■), the Poisson 
bracket is dehned by 




£(s)'-(s;a- 


( 2 . 11 ) 


where “t” denotes the matrix transpose. This bracket obeys the standard antisymmetry 
and Jacobi properties [45], due to the skew-adjoint form of P combined with the property 
that the components of F have no dependence on the variables v and w. One immediate 
consequence of a Hamiltonian formulation is that the Hamiltonian is a constant of motion, 
namely 


m 

dt 


= 0 


( 2 . 12 ) 


whenever v and w obey suitable asymptotic decay conditions for large x. 
The Hamiltonian for the highly nonlinear wave equation fll.ip is given by 






)■ 


(2.13) 


This Hamiltonian will physically describe a conserved energy, consisting of a kinetic energy 
term plus a potential energy term. 

It may also be of interest to note that this Hamiltonian formulation can be re-expressed as a 
second-order equation in time derivatives for the strain v. Consider the general Hamiltonian 
system fl2.8p - fl2.10p . By taking the time derivative of the v equation and substituting the w 
equation, we obtain Vu = —Wtx = D‘^{^'H/6v). This yields a wave equation with the form 


{^Hq/6v'^xXJ Hq Lqj 


(2.14) 


which is neither Lagrangian nor Hamiltonian in terms of v. Here Flo can be interpreted as 
the potential energy. For the highly nonlinear wave equation fll.ip . we have 

= + (2.15) 

which gives 

+ lR‘'v'"~"vl)xx- (2.16) 

We will next derive the conservation laws and corresponding conserved quantities for the 
highly nonlinear wave equation fll.ip . 


2.2. Conservation laws. The most commonly used method for hnding conservation laws 
of Lagrangian differential equations is Noether’s theorem [4^146] . This method produces a 
conservation law from every symmetry of a Lagrangian. For example, if a Lagrangian is 
invariant to within a total divergence under time translations, then this symmetry typically 
produces the conservation law for energy. Some drawbacks to using Noether’s theorem are 
that, hrst, the symmetries of the Lagrangian need to be determined, and second, the variation 
of the Lagrangian under each symmetry needs to be found to obtain the expression for the 
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resulting conservation law. We will instead use an alternative, modern version [161I17] of 
Noether’s theorem, employing multipliers rather than symmetries, which avoids the usual 
drawbacks and which is simpler to carry out. 

To explain the method and to review the theory of conservation laws, we consider a general 
fourth-order wave equation Uu = F{u,Ux,Uxx,Uxxx,Uxxxx)- A conservation law for a wave 
equation of this form consists of a local continuity equation 

DtT + DxX = 0 (2.17) 

holding for all solutions u{t,x) of the wave equation, where T is the conserved density and 
X is the spatial flux, which are given by functions of t, x, u, Ut, and x-derivatives of u and 
Ut- Note that Uu (and any derivatives of it) can be assumed to be eliminated from T and X 
through the wave equation. 

Every conservation law fl2.17p yields a conserved quantity dehned by the global continuity 
equation 

^[Tdx = -X (2.18) 

dt an 

on a given spatial domain with endpoints dfl. For solutions u{t,x) of the wave equation 
having suitable boundary conditions at the endpoints of the domain, the spatial flux terms 
will vanish so then the quantity 


C[u\ 


T dx = const. 


(2.19) 


Jn 

is a constant of motion. Two conservation laws are physically equivalent if they give the 
same conserved quantity up to boundary terms. This occurs if (and only if) their conserved 
densities T differ by a total x-derivative, ■, and correspondingly, their fluxes X differ by 
a total time-derivative, for all solutions u{t,x) of the wave equation. Consequently, 

a conservation law will be physically trivial if, on all solutions, the density and flux have 
the form T = and X = —Dt"^ whereby the local continuity equation (12.171) will hold 
trivially. Thus, equivalent conservation laws differ by a trivial conservation law. 

Any conservation law fl2.17l) has an equivalent characteristic form that holds as an identity 
obtained by moving off solutions of the wave equation: 


DtT + Dx{X-<!>) = {utt-F)Q, 


( 2 . 20 ) 


where 

^ = {utt — F)Eut^(T) + Dx{utt — F)Eut^^{T) + ■ ■ ■ ( 2 . 21 ) 

is a trivial flux which vanishes on all solutions of the wave equation, and where 


Q = EuAT) 


( 2 . 22 ) 


is called the multiplier. Here 


F = — -n — 


d: 


d 


dZx 




(2.23) 


denotes the (spatial) Euler-Lagrange operator with respect to a variable Z. This operator 
has the property that Ez{f) vanishes precisely if / is a total x-derivative of an expression 
in terms of x, Z, and x-derivatives of Z. We provide a short derivation of equations (I2.2np - 
(I2.22P at the end of this section. 
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Note that the characteristic equation fl2.20l) reduces to the conservation law fl2.17p when 
u satishes the wave equation uu = F. This formulation of conservation laws is important for 
several reasons. First, it shows that there do not exist any non-trivial conserved densities 
T that depend only on t, x, u, and a:-derivatives of u, Ut- In particular, T must have some 
essential dependence on Ut (or its a:-derivatives), as otherwise D^T + D^X cannot vanish 
whenever u satishes Uu = F. Second, from this conclusion, it follows that if Q = Eu^(T) = 0 
then T = T + holds where T has no dependence on ut (and its x-derivatives), thereby 
implying T = 0. This shows that the multiplier of a conservation law will be trivial precisely 
when the conserved density is trivial. Thus, there is a one-to-one relationship between 
non-trivial conserved densities (up to equivalence) and non-vanishing multipliers. Moreover, 
when Q ^ 0, the conserved density T and the hux X can be recovered from Q by a direct 
integration of the characteristic equation fl2.20p . 

This reduces the problem of hnding conservation laws to the simpler problem of hnding 
multipliers. From the characteristic equation fl2.20p . we see that a function Q will be a 
multiplier if (and only if) its product with the wave equation has the form of a total t- 
derivative plus a total x-derivative. Such expressions have the characterization that their 
variational derivative with respect to u vanishes identically. This condition 

r 

— ((Uii-F)Q) =0 (2.24) 

will split into a system of linear partial differential equations that determine Q. First, from 
relation fl2.22p . all multipliers are functions only of f, x, u, Ut, and x-derivatives of u and Ut- 
Then, since Uu, Uttx, etc. do not appear in Q, the coefficients of these variables in the condi¬ 
tion fl2.24|) must vanish separately. The resulting equations give an overdetermined system 
that can be straightforwardly solved for Q. We emphasize that this system is applicable 
whether or not a Lagrangian formulation exists for the wave equation Uu = F. 

In the case when a wave equation uu = F{u,Ux,Uxx,Uxxx,Uxxxx) possesses a Lagrangian 
formulation fl2.3p . the determining system for multipliers Q can be simplihed using some 
tools from variational calculus. This leads to the system 


DfQ = F^Q + F„,D,Q + F„,,DlQ + F^...DlQ + Fu..„DiQ, (2.25) 

Q^.=E^,{Q), =-£■«((?), Q„,„ = BS’(Q). .... ( 2 . 26 ) 


where 


A = A 


d d ^ d d 


utt=F 


n 17 ^ ^ 

F)xF ^ 'i^txx " 


(2.27) 


' dux dutx duxx 

is the total time derivative evaluated on the space of solutions of the wave equation, and 
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+ • • • , 


(2.28) 


are higher-order Euler-Lagrange operators usiiiii. The hrst equation in the system fl2.2 
(I2.26P coincides with the determining equation for inhnitesimal symmetries X = Qd/du of the 
wave equation, with the symmetry dehned to act on solutions u{t,x) of the wave equation 
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by the infinitesimal transformation u —?• u + eQ + 0{e^) in terms of a parameter e. The 
remaining equations in the system fl2.25p - fl2.26p comprise the Helmholtz conditions [15] which 
are necessary and sufficient for Q to have the form of an Euler-Lagrange expression fl2.22p . 
The Helmholtz conditions turn out to be equivalent to the condition that the Lagrangian of 
the wave equation is invariant to within a total derivative under the prolonged symmetry 
prX = Qd/du + D^Qd/du^ + • • •. Consequently, multipliers correspond directly 

to symmetries of the Lagrangian. This provides a simple formulation of Noether’s theorem 
that avoids any use of the Lagrangian itself [^6] - H8] . 

For the highly nonlinear wave equation fll.ip . we will now present all conservation laws 
given by conserved densities of the first-order form 

T{t,x,u,u^,u^^,ut) (2.29) 

which can be expected to encompass all of the physically relevant conserved quantities m- 
The corresponding multipliers will have the form 


Q{t,x,u,u^,ut) (2.30) 

as shown by the variational relation fl2.22l) . The results are obtained by directly solving 
the determining system fl2.25l) - fl2.26l) . For this class of multipliers fl2.30l) . note that the 
characteristic equation simplifies to 


DtT + D,X = {uu-F)Q 


(2.31) 


since the trivial flux term fl2.21l) vanishes. We use this equation to obtain the conserved 
density and spatial flux corresponding to each multiplier. Specifically, as Uu, Uuxi etc. do 
not appear in T, X, and Q, the coefficients of these variables in equation fl2.3ip must vanish 
separately, whereby this equation splits into a system of linear partial differential equations 
that are easily integrated to find explicit expressions for T and X. 

All calculations for our results have been carried out using Maple. 

The first-order multipliers admitted by the highly nonlinear wave equation fll.ip are given 
by 

Qi = 1, Q2 = t, Qs = -Ux, Qi = Ut. (2.32) 

The corresponding conserved density and flux determined by each multiplier are given by, 
respectively. 




X Kx + iK)-^ 
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As = \ul + c^{\R^uluxxx + 
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-2 T^2~k-l~2 




(2.33) 

(2.34) 

(2.35) 
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These conservation laws yield the following conserved qnantities: 

Ut dx, 


Ci[u\ = 
C2[u\ = 
Cslu] = 


tut — udx = tCi [u] — u dx, 
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—UtUx dx, 

Jn 

CM = [ {¥ 
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)) dx. 


(2.37) 

(2.38) 

(2.39) 

(2.40) 


The qnantity (I2.4np coincides with the Hamiltonian (I2.13p which describes the physical 
total energy of wave solntions u{t,x), while the qnantity (I2.39p describes the physical total 
momentnm of wave solntions. These interpretations are reinforced by the observation that 
the respective mnltipliers correspond to an inhnitesimal time-translation X = Utd/du and 
space-translation X = Uxd/du from the Lagrangian viewpoint. 

The qnantities fl2.37p and fl2.38p have a less familiar physical meaning. We hrst note 
u dx can be viewed as the mean amplitnde associated to a wave solntion u{t, x). Then we 
see 


f ~ j dCi[u] 

— udx = —-— = 0 
dt^ L dt 


(2.41) 


holds for wave solntions having snfficient asymptotic decay for large x. Thus the mean 
amplitude u dx obeys the equation of free particle motion. Moreover, at the initial time 
t = 0, the mean amplitude is related to the quantities fl2.37p and fl2.38p by 

= ‘‘ 


u dx 
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= Mu] 


(2.42) 


Hence these quantities are just the initial value of the mean amplitude and the initial value 
of the mean longitudinal-momentum. 

2.3. Derivation of multiplier equation. Here we give a short derivation of the char¬ 
acteristic equation fl2.20p for conservation laws of a general fourth-order wave equation 
uu = F{u, '^XXI '^XXXj "^XXXX^' 

Since any conservation law fl2.17p holds for all solutions of Uu = F, we can re-write it 
using the time derivative fl2.27p evaluated on solutions: 


Df T DxX — 0 . 


(2.43) 

Then, as Uu does not appear in this equation, it holds with u{t, x) replaced by an arbitrary 
function. We now substitute the relation 


o o O 

Dt-Dt = {F- uu)j— + Dx{F - uu)j— + Dl{F - uu)j: - 

OUt UUtx UUtxx 


+ ••• 


which yields an identity 

DtT + DxX = {uu - F)Ko + Dx{uu - F)K, + Dl{uu - F)K2 + 
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Next we use integration by parts on each term on the righthand side in the identity fl2.45p : 
D^{utt - F)Ki = -{utt - F)D^Ki + D^[{utt - F)Ki), 

Dliutt - F)K2 = {utt - F)DIK2 + D^{D^{uu - F)K2 - {uu - F)D^K2), (2.47) 


Collecting the total derivative terms, we get 

DtT + D^X = {utt - F){Ko - + DIK 2 + ■ ■ ■) 

+ D^{{uu - F){K, - D,K2 + ■■■) + -F){K2 + ---))+--- . 

(2.48) 

After substituting expressions fl2.46p into this equation, we obtain the characteristic equation 
fl2.20p together with the expressions 02.221) and 02.2ip for the multiplier and the trivial flux 
term. 


3. Travelling wave solutions 


For a general fourth-order wave equation uu = F{u,Ux,Uxx,Uxxx,Uxxxx), the form of a 
traveling wave is u = /(^) with ^ = x — ut, where the constant z/ is the velocity of the wave, 
and the function / is the wave prohle. Since Ux = f and Ut = —I'f = —i^Ux, where a prime 
denotes the derivative with respect to the wave equation will reduce to a fourth-order 
ODE 

= F{f,f\f\r\r) ( 3 . 1 ) 

for /(O- 

We will now explain how to use conservation laws to directly provide hrst integrals of the 
traveling wave ODE and also to obtain conserved quantities for the solitary wave solutions 
and nonlinear periodic solutions of this ODE. 

First, suppose a conservation law DtT + DxX = 0 of the wave equation Uu = F does not 
contain the variables t and x explicitly. Then this conservation law will give rise to a hrst 
integral of the traveling wave ODE fl3.1l) by the reductions 


-=/(«) I-=/(«) de’ 

yielding 

4((x-.t)L^,„)=o. 

The hrst integral is thus given by 

X — vT = C = const.. 


(3.2) 


(3.3) 

(3.4) 


It has the physical meaning of the spatial hux in a reference frame moving with speed z/ 
(namely, the rest frame of the travelling wave). Existence of three functionally-independent 
hrst integrals fl3.4|) will then allow the traveling wave ODE fl3.ip to be reduced to a hrst-order 
ODE, f = I{f). 

Second, suppose a conservation law DtT + DxX = 0 has vanishing hux at the endpoints of 
a spatial domain D in a: for a travelling wave solution u = f{^) of the wave equation Uu = F. 
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(3.5) 


Then the conservation law yields a conserved quantity 

^ f T\ ^ I T\ , = Q 

dtJn dtJn 

where the domain will be = (—C) 0 , oo) for solitary wave solutions, or ^2 = (—L/2, L/2) for 
nonlinear periodic wave solutions with wavelength L > 0. If the conserved density T does 
not contain the variables t and x explicitly, then will be an expression involving only 

/ and its derivatives up to fourth order (i.e. ^ will not appear explicitly). We can simplify 
this expression by using the reduced ODE /' = /(/) to eliminate all derivatives of /, so that 
^U=/(f) = T(/). If the wave prohle of a travelling wave solution u = f{^) on the domain 
D has a single maximum with an amplitude /max at the center .^ = 0 of the domain and 
a minimum with an amplitude /min at the domain endpoints ^ = dfl, then the conserved 
quantity for this solution can be expressed as 

assuming that the wave prohle is symmetric around / = 0. This expression can be evaluated 
without the need for the explicit form of the solution /(/). 

3.1. Highly nonlinear wave eqnation. We will now apply the method outlined in sec¬ 
tion [ 1 ] to derive all travelling wave solutions for the highly nonlinear wave equation fll.ll) . 
To begin, we substitute the traveling wave expression 

u = /(/), ^ = X — ut, 1 / = const. (3.7) 

into the wave equation fll.ip . This yields the fourth-order ODE 

(^/c) + «(/')"■"(/'? + (3.8) 

where c, a, [3, 7 are given by expressions fll.2l) - fll.3p . and where the nonlinearity exponent 
k obeys [TTfH3] 

k>l. (3.9) 

The physical variable that will physically describe a propagating wave is the strain v = 
Ux = —Ux > 0 (which is positive because the underlying discrete dynamical system has no 
restoring forces). Hence, we require the condition 

/' > 0 (3.10) 

on the travelling wave prohle /(/). Moreover, since only derivatives of f{^) appear in the 
fourth-order ODE fl3.8l) . we will rewrite it in terms of the strain variable 

v = ~g = f>0 (3.11) 

which gives a third-order ODE 

{u/crg' = f-^~g' + at-\g'f + • (3.12) 

We now seek to hnd all solitary wave solutions and all periodic wave solutions of the strain 
ODE flTT^ . 

The hrst step in our method is to make use of the conservation laws fl2.33p - fl2.36p derived 
in section [2] for the highly nonlinear wave equation fll.ll) . By applying the direct reduction 
approach fl3.4l) . we obtain hrst integrals for the ODE fl3.12l) from the conservation laws fl2.33p . 
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fl2.35p . and fl2.36p . These three conservation laws, which do not contain t and x explicitly, 
respectively describe the mean longitndinal-momentnm, the total momentnm, and the total 
energy for solntions of the wave eqnation fll.ip . 

The first integral arising from the mean longitndinal-momentnm fl2.33p is given by 

Cl = {u/cf~g - If - = const.. (3.13) 

Physically, it measnres the spatial flux of mean longitudinal-momentum in the rest frame of 
the travelling wave. Both the total momentum (I2.35p and the total energy (12.361) yield the 
same first integral (up to a multiplicative constant) 

C 2 = - ^R^f-^~g'^ - \R^f~g'' = const.. (3.14) 

Physically, this hrst integral measures the spatial flux of total energy/momentum in the rest 
frame of the travelling wave. We use the hrst integral fl3.13p to solve for f in terms of g' 
and g, and substitute this expression into the hrst integral fl3.14p . This immediately yields 
the hrst-order separable ODE 

{g'f = - Ci~g + C2- (3.15) 

We write this ODE in dimensionless form by introducing the scaling 

g = ^g, i = (3.16) 

Cl =-CiA;(fc + l)A-^ C 2 = C 2 k{k + (3.17) 

If we choose 

X= {fu/cyk{k + l)Y^^''~^\ u= ^J\k{k + l)R, (3.18) 

then the ODE becomes 

{g'f = G{g), (3.19) 

with 

G{g) = g^ ^{g"^ + Cig + C2 — g^^f, g>^ (3.20) 

where ^ and g{f are now dimensionless, and where Ci and C 2 are dimensionless arbitrary 
constants. 

The next step in our method is to classify all solutions of the hrst-order separable ODE 
fl3.19p by employing a modihed energy analysis. We hrst write C 2 = C, where E represents 
a dimensionless energy for a solution g{C). We view Ci as a free parameter on which the 
solutions g{f will depend. Then we introduce the potential energy function 

V{g)=g^^^-g{g + C^), (3.21) 

and write the ODE in the form 

gk-ig'2 + ^ E ( 3 . 22 ) 

which is a sum of a kinetic-type energy term and a potential energy term. The main part of 
the analysis is that we need to understand the shape of the potential fl3.2ip as a function of 
g > ^. To begin, we note its hrst and second derivatives are given by: 

V'{g) = (1 + k)g^ -2g- Ci, V"{g) = kf + k)g^-^ - 2. (3.23) 


We then see that there is only one point at which the convexity of V{g) changes: 


g 


* 


/ 2 0/(*:-l) 

\k{k+l)) 
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(3.24) 

























Since k belongs to the range fl3.9p . note that we have 


0 < / < 1. (3.25) 

We now hnd all of the critical points g = of V{g). The number and nature of these 
points will depend on the value of Ci. First, we observe that hnding the roots of V'{g) = 0 
is equivalent to hnding the intersection points of the curve (1 + k)g^ and the line 2g + Ci. 
Now consider the slope of the curve, which is given by (1 + k)kg^~^. At 5 ^ = 0, the slope 
vanishes. For > 0, the slope is increasing and becomes equal to the slope of the line at 
g = g*. For g > g*, the curve has a greater slope than the slope of the line. Hence, if 
Cl > 0, then for g > 0 the curve will intersect the line at exactly one point g = gl, which 
will he to the right of g*. Since V"{gl) > 0, we conclude that gl is a local minimum of 
V{g). Next, if Ci < 0, then there are three cases, depending on the sign of the difference 
(1 + k)g*’‘ - {2g* + Ci) = - Ci where 

Cl = {l + k)g*^-2g* =<f). (3.26) 

If the sign of Cl — C\ is negative, then for g > 0 the line will intersect the curve in two 
distinct points g = g\ and g = g^ such that 0 < gl < g* < g^. The critical point at gl is 
a local maximum of V{g), while the other critical point g^ is a local minimum of V{g). If 
Cl — Cl is zero, then for g > 0 the line will intersect the curve only at the point (i.e., as a 
tangent line) g = g*, which is an inflection point of V{g). If the sign of Cl — Ci is positive, 
then for g > 0 the line will not intersect the curve. This completes the determination of all 
of the critical points of V{g). 

Finally, we look at the asymptotic behaviour of V{g) as g ^ 0 and g —?• cxd. For large g, the 
dominant term in V{g) is g^^^ 3> g{g + Ci), so then V{g) —>■ 00. For g near 0, the dominant 
term in V{g) is either \Ci\g S> g^~^^ if Ci 7 ^ 0, or g"^ S> g^^^ if Ci = 0. Combining this 
asymptotic analysis with the critical point analysis, we obtain a full picture of the potential 
energy function V{g). In particular, there are hve distinct cases, depending on the value of 
Cl, as illustrated in Fig. [1] 

We can now classify the different types of solutions of ODE fl3.22p by using an energy 
analysis based on the features of the energy potential V{g). Consider the intersection of 
an energy line E = const, and the potential V{g) for g > f). The ODE fl3.22p shows that 
gk-igi 2 _ g intersection point, and this implies g' = 0 (since g > 0). Hence, 

all such intersection points are critical points of the wave prohle g{^)- This leads to four 
different types of solutions, two of which have the feature that g{^) vanishes at isolated 
points 

Additional types of solutions arise if we also consider intersection points given by = 0, 
which can occur (only) if E = 0. Compared to the four main types, these additional solutions 
exhibit different behaviour at the points ^ where g{^) = 0 . 

A summary of all of the different solution types is stated in Table [T] Details will be derived 
in subsequent subsections 13.2113.61 

The presence of a node in a travelling wave g{^) means that the strain fl3.1ip is no longer 
positive but vanishes at isolated points ^o- This can be viewed as a limiting case of vanishingly 
small strain. 
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Figure 1. Potential energy function V{g) = — g{g + Ci) for the Hertz 

potential k = 3/2. 


To understand how the different types of solution arise, we integrate the hrst-order sepa¬ 
rable ODE (13.221) to get 


± 



(k-l)/2 

. = dg = ^ 

VE-V{g) 


(3.27) 


where (70 > 0 is a root of V(g) 


= E. Then we can factorize the denominator expression 
E - V{g) = {g- goTA^g) 
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(3.28) 




























Table 1. Types of travelling wave solutions 


Intersection points 
V{g) = E 

Type of point 

9 

Nonlinearity 

k 

Solution type 

9{.i) 

Domain 

two g > 0 

one local maximum 
one non-critical 

> 1 

solitary wave 

(—oo, oo) 

one g > 0 

local maximum 

> 1 

cusp-nodal 
solitary wave 

(—oo, cxd) 

two g > 0 

non-critical 

> 1 

periodic wave 

—L T7\ 

. 2 ’ 2 J 

one g > 0 

non-critical 

> 1 

cusp-nodal 
periodic wave 

— -1 

1. 2 ’ 2 J 

one g = 0 
one g > 1 

non-critical 

non-critical 

< 2 

nodal periodic wave 

-1 

one g = 0 
one g > 1 

non-critical 

non-critical 

= 2 

corner-no dal 
periodic wave 

L 2 ’ 2 J 

one g = 0 
one g > 1 

non-critical 

non-critical 

> 2 

cusp-nodal 
periodic wave 

k] 

y 21 2i 

one g = 0 
one g = 1 

critical 

non-critical 

< 3 

nodal periodic wave 

y 21 2i 

one g = 0 
one g = 1 

critical 

non-critical 

= 3 

corner-nodal 
periodic wave 

k] 

1 . 2 ’ 2 J 

one g = 0 
one g = 1 

critical 

non-critical 

> 3 

cusp-nodal 
periodic wave 

L 2 ’ 2 J 


where m is some positive integer (i.e., the degeneracy of the root) and A{g) is some function 
of g such that A{go) is non-singular. This factorization will hold under the conditions that 
E — V{g) and its derivatives up to order m — 1 vanish at g = go. Note, from the shape of 
V{g), the maximum possible degeneracy is m = 3, and the root go is a critical point of V{g) 
iff m 7 ^ 1. Next we expand the integral fl3.27p asymptotically for g near go > 0: 
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{k-l)/2 


VE - V{g) 


dg 


96 


(fc-l)/2 


dg 


VlMMJ \9-9o\^/^ 

( fc - i )/2 I 2 | 5 ' — 


96 


V\M9o)\ 



- 1/2 


m = 1 
m = 2 
m = 3. 


(3.29) 


In the case m = 1, this integral converges, and so g{^) will equal go at a hnite value ^ = .^o- 
In the other two cases, the integral diverges, whereby g{^) will asymptotically approach go 
as |.^| gets large. This asymptotic tail will have exponential decay of \g — go\ —t 0 in the 
case m = 2 and power decay oi \g — go\ —)■ 0 in the case m = 3. The global behaviour of 
the solution g{^) as a function of ^ depends on whether the domain for the integral fl3.27p is 
given hy g > go OI g < go, as determined by the condition V{g) < E so that ^yE — V{g) is 
well-defined. 

When V{g) < E holds for some g > go > the shape of V{g) shows that Ci must be 
greater than Cj, and that V{g) = E has a second root gi. This root satishes gi > go > 0 
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and is not a critical point of V{g). Conseqnently, the domain for the solntion integral fl3.27p 
will he go < g < Qi, and the resnlting solntion g{^) will have the properties that it will reach 
gi at a hnite valne of and that it will reach go either at a different hnite valne of ^ if go is 
not a critical point of V{g), see Fig. [2] or asymptotically for large |^| if go is a critical point 
of V{g), see Fig. [31 Moreover, in the latter case, the shape of V{g) shows that V"{go) ^ 0, 
so the degeneracy of go will be m = 2. These two possibilities for the global behavionr of 
g{^) thereby describe, respectively, a periodic wave (obtained by piecing together translated 
copies of the wave prohle) whose maximnm and mininnm amplitndes are gi and go, and a 
solitary wave with a maximnm amplitnde gi and an exponential tail decaying to go- 


V(g) V(g) 




Figure 2. Conditions for a periodic wave solntion for the potential energy 
fnnction V{g) = — g{g + Ci) for the Hertz potential k = 3/2. 


V(g) 



Figure 3. Conditions for a solitary wave solntion for the potential energy 
fnnction V{g) = — g{g + Ci) for the Hertz potential k = 3/2. 


When V{g) < E holds for some g < go > 0 with E > 0, the shape of V{g) shows that 
V{g) = E > 0 has no roots other than go and that there is no restriction on Ci- Conseqnently, 
the domain for the solntion integral (I3.27p will he 0 < g < go- The resnlting solntion g{^) 
will reach 0 at a valne of ^ determined by the convergence of the integral (13.271) near g = 0. 
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dg = 0 


(3.30) 


In particular, due to fc > 1 and E > 0, we have 

,(fc-l )/2 


9 ^ 


VE - V{g) 


dg 


^/E 


9 


(fc-l )/2 


and 

{g'r ~ Eg^-’^ = oo. (3.31) 

Hence, g{^) will go to 0 at a finite value of while g'{^) will blow up. This describes a 
nonlinear wave g{^) that has a nodal cusp. Furthermore, this solution will reach go either 
at a hnite value of ^ if go is a not critical point of V{g), see Fig. 01 or asymptotically for 
large |.^| if go is a critical point of V{g), see Fig. [5l These respective cases describe a cusped 
periodic wave (obtained by piecing together translated copies of the wave prohle), and a 
cusped solitary wave. In the latter case, the shape of V{g) shows that either V''{go) 7 ^ 0 if 
Cl > C^, whereby the degeneracy of go will be m = 2 , or V''{go) = 0 if Ci = C^, whereby the 
degeneracy of go will be m = 3. These two possibilities correspond to the tail of the cusped 
solitary wave having either exponential decay or power decay of — ^fol —t 0 for large |.^|. 



V(g) 




E>0 


E>0 

g 


Figure 4. Conditions for a cusped periodic wave solution for the potential 
energy function V{g) = g^~^^ — g{g + Ci) for the Hertz potential k = 3/2. 


A limiting case occurs when V{g) < E = {] holds for some g > go = where H(0) = 0 
and H'(0) = —Ci. In this case, as shown by the shape of V{g), see Fig. 01 we have Ci > 0, 
and V{g) = C = 0 has a second root gi > g* > Q which is not a critical point of V{g). 
Moreover, we see 5^0 = 0 is a critical point of V(g) iff Ci = 0. The domain for the solution 
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V(g) V(g) 




Figure 5. Conditions for a cusped solitary wave solution for the potential 
energy function V{g) = — g{g + Ci) for the Hertz potential k = 3/2. 


integral fl3.27l) thereby will he 0 < g < gi, and the resulting solution g{^) will reach gi > 0 
at a hnite value of The solution g{^) will reach = 0 at a value of ^ determined by 
expanding the solution integral (13.2711 asymptotically for g near 0: 

M g{k-l)/2 ra ^^2-1 J^^(fc-i)/2 _ 0, ^1=0 

Jo VE-V{g) Jo ( 3 - 32 ) 

~ 0 


which is due to k > 1. Thus, since this integral converges to yield 0, we can express the 
solution g{^) as 


rg gk/2-l 

Jo ^/Ci + g- g’^ 

where ^'(O) = 0. The ODE (13.221) then shows that 


dg 




(3.33) 


This implies 


and 


id') 


)2 ^ 1 

g3 

C,g^-\ 

Ci = 0 

as 

Cl > 0 

^ ^ 0. 

(3.34) 

fo. 

k <3 




r’ 

k = 3 

when F'(0) = 

= 0, Fig.EKi), 

(3.35) 

[oo. 

k > 3 




0 , 

k<2 




F'(0), 

k = 2 

when F'(0) 

o' 

V 

(3.36) 

oo. 

k>2 





Therefore, the solution g{^) describes a periodic wave in which = 0 modulo L is a node that 
exhibits a cusp, a corner, or a local minimum, depending on the values of k and F'(0) = —Ci, 
where L is the wavelength. (Note that the case F'(0) = 0 distinguishes when the set of 
constant strain solutions, g = const., includes g = 0.) 
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V(g) 


V(g) 




Figure 6 . Conditions for the limiting case F = 0 (nodal, corner, cusp) peri¬ 
odic wave solutions for the potential energy function V{g) = — g{g -|- Ci) 

for the Hertz potential k = 3/2. 

We remark that the presence of a node in a travelling wave g{^) means that the strain 
fl3.1ip is no longer positive but vanishes at isolated points ^o- This can be viewed as a limiting 
case of vanishingly small strain. We also remark that the wave profile g{^) for cusped waves 
and peaked waves is not a differentiable function of ^ at the locations of the cusps and peaks. 
Thus, at these values of g{^) does not satisfy the ODE fl3.22p . As a consequence, for g{^) 
to make sense as an actual solution, we would need to seek a suitable weak formulation of 
the ODE, similarly to what is done for peakon wave solutions |STHS3] of the Camassa-Holm 
equation and other related nonlinear wave equations. 

In the next step of our method, based on the preceding analysis of the solution integral 
fl3.27p . we determine the conditions on the constants E, go, and Ci = F'(0) under which 
the various types of solutions g(^) (listed in Tabled]) are produced. For carrying out this 
analysis, it will be useful to introduce the notation 

— If- 

Sn{a,b) = Sn{b,a) = -—, a ^ b, n>0. (3.37) 

a — b 

As a function of a and b, Sn{a,b) has the following properties that we will need. First, for 
a > b > 0 and 6 > a > 0, S'„(a, b) is non-singular when n > 0; daSn{a, b) and d^Sn^a, b) are 
non-singular when n > 1 and n > 2 , respectively. Second, Sn{a,b) obeys the inequalities 


Sn{a,b) <n , , , n 

a + b 

IV 

o 

(3.38) 

and 

Sn{a,b) > Sn{a,c) > 0, 

n > 0 

(3.39) 


n > 1 

(3.40) 


n > 1 

(3.41) 

dlSn{a,b) > dlSn{a,c) > 0 , 

n > 2 

(3.42) 

for a > 6 > c > 0. Third, S'„(a, b) has the limits 

•) ^^\b=a 

|n(n — l)a"' ^ 

(3.43) 
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for a > 0 . 


(3.44) 


3.2. Periodic waves. To obtain a periodic nonlinear wave solntion g{^) > 0, we must have 

(Ti > 0 : 0> E> 

0 > Cl > C* : V{gt) >E> Vig^,) 

E = V{go) = V(gi), 0 < go < gi (3.45) 

V'(go) < 0 , V'(gi) > 0 (3.46) 

V{g) <V(go) for go < g < gi (3.47) 


where V'{gl) = V'{g 2 ) = 0. 

The second condition fl3.45l) yields E = gQ~^^ — go{go + Ci) = g^^^ — gi{gi + Ci). These 
two equations determine 


Cl — Sk+i{gi, go) — {gi + go) 

E = 5'i5'o(l - Sk{gi,go)) 

in terms of go and gi. Then we have 

y{g) = gig'" - g + gi +go- Sk+i{gi,go)) 

which gives 

E-V{g) = {gi - g){g- go)A{g) 

where 


A{g) = 


Sk+i{g,gi) - Sk+i{g,go) 


^ ^ giSk{g,gi) - goSk{g,go) _ ^ 
gi — go 


(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 


gi — go 

Next, we hnd that the third condition fl3.46p yields 

gi{Sk{gi,go) - i) > go{kgo~^ - 1) (3.53) 

and 

gi{kgi~^ - 1) > go{Sk{gi,go) - i)- (3.54) 

Comparing these two inequalities, we see that the second one follows from the hrst one, 
because k{g'l + g^) - (gi + go)'Sk(gi, go) > 0 due to property (13.381) . using 

0 < < gi- (3.55) 


Then, the remaining conditions fl3.47p and fl3.44p will hold because V(go) — V(g) is increasing 
for both g — go and gi — g increasing up to the single critical point g = g^ between go and 
gi. Finally, we note that inequality fl3.53p imposes a simple lower bound on gi as follows. 
Consider the extended inequality gi{Sk{gi-i go) — 1) > go{kgo~^ — 1) with Q < go < gi- 
This inequality can be re-written as Bo{go) + Bi{go)gi > C{gi), which says that the line 
Bo{qo) + - 81 ( 90)91 lies above the curve C( 9 i), where Boia) = q‘^{kq^~^ — 1 ), BAq) = q{2 — 
{k + l)g’^-^), C{g) = g‘^{l-g^~^). We see that ii go = gi then Bo{gi) + Bi{gi)gi = C{gi) and 
B\{gi) = C'{gi), so thus the line touches the curve as a tangent line. Hence, the boundary of 
inequality (13.531) consists of the points 5^1 > 0 at which the line Bo{go) + Bi{go)gi is tangent 
to the curve C{gi). The minimum gi in this set of tangent points can be seen geometrically 
to come from the critical point ^fg of the set of lines Bo{go) + Bi{go)gi as given by BQ{gQ) = 0 
and B[{gQ) = 0. We have BQ{g) = k{k -|- l)g^ — 2g and B[{g) = 2 — k{k + l)g^~^, which 
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yields = g*, where g* is expression fl3.24p . This gives gi = g* for the minimum of gi in 
the set of tangent points. We thus have the lower bound 

91 > 9 * (3.56) 

under which the inequality fl3.53p will hold. 

Therefore, a periodic wave exists iff the constants go and gi satisfy the inequalities fl3.53p . 
fl3.55p . fl3.56p . Then the solution integral fl3.27p is given by 

rg a(k-l)/2 

± / /, „ ofT ‘‘3 = (=f i/2 (3.57) 

“'90 y [91 — 9){9 — 9o)A{9) 

with 9 o ^ 9 ^ 9 i, where 

f9l g{k-l)/2 

L = 2 , ^ dg>0 (3.58) 

Jgo Vi9i- 9) {9 - 9o)M9) 

is the wavelength. The wave prohle g{^) for this solution fl3.57p will reach a local maximum 
at .^ = 0 modulo L and a local minimum go a.t ^ = ±L/2 modulo L. 


3.3. Cusped periodic waves. This type of periodic solution g{^) > 0 arises when 



U > 0 

(3.59) 


E = V{go) 

(3.60) 


V'{go) > 0 

(3.61) 


V{g) < V{go) foT0< g < go- 

(3.62) 

The second condition fl3.60p vields E = — aAon + Uih which determines 



Cl — go — 9o — E/go- 

(3.63) 

This gives 

V{ 9 ) = { 9 ^ ~ 9o ~ 9 + 9o + E/go)9 

(3.64) 

and thus we have 

E -V{g) = {go- 9)A{g) 

(3.65) 

where 




= Sk{g, go)g -g + E/go = Sk+i{g, go)- g- go+ E/go- 

(3.66) 

From the hrst and third and conditions fl3.59p and fl3.61p. we have 



E/go > 0, 

(3.67) 


E/go > go{l - kg^-^). 

(3.68) 


The fourth condition fl3.62l) reduces to A{g) > 0. We note ^4(0) = E/go > 0 and A{go) = 
kgo — go + E/go > 0 hy using inequalities fl3.67l) - fl3.68l) and property fl3.43p . Now consider 
A'{g) = dgSk+i{g, go) — 1 from expression fl3.66l) . Using properties fl3.40l) and fl3.43p . we 
see A'{g) < A'(go) = dgSk+i{go, go) - 1 = {go/ 9 *)''~^ - ^ ior 0 < g < go, from expressions 
fl3.66p and fl3.24p . Hence, we obtain A'{g) < 0 if (70 < 9 *: which then implies A{g) > A{go) 
ioY 0 < g < go- Since ^(^fo) > 0, we conclude A{g) > 0 for 0 < (7 < (70 in fhe case 
9 o ^ 9 *■ The complementary case go > g* is more complicated and leads to two subcases as 
follows. We first note A'{g) > dgSk+i{g, g) — 1 = |(fc + l)kg’^~^ — 1 = {g/g*)’^~^ — 1 from 
properties fl3.4Up and fl3.43l) . This implies ^'(^f) > 0 for 77 > whereby ^(^f) is increasing 
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for g > g*. Also, we have A"{g) = dgSk+i{g, go) > 0 by property fl3.42p . Next we note 
that A'(0) = dgSk+i{g, go)\g=o — 1 = gQ~^ — 1 by direct evaluation. This gives a case split, 
depending on whether g^ > 1 or g^ < 1. If g^ > 1, we have A'(0) > 0, which then implies 
that A{g) is increasing for 0 < (/ < go, since A''{g) > 0. We thereby conclude that ^(^f) > 0 
due to A(0) > 0. If instead go < 1, then we have A'(0) < 0, which implies A{g) has a unique 
minimum for some g = g^ < S'*, since A'{g) > 0 for g > g* and A"{g) > 0 for 0 < g < go- 
This minimum is determined by 

dgSk+i{g,go)\g=g* = 1. (3.69) 

Hence, we conclude A{g) > Ai^g^), and then we require A^g^) > 0 so that A{g) > 0 will hold 
for 0 < S' < S'o- This yields the condition E/gQ > (1 — *S'fc(s'o,S'o))S'o from expression fl3.66p . 
Finally, we note that equation fl3.69p can be expressed in an equivalent form by expanding 
out the derivative, which gives *S'fc+i(s'o, S'o) = (k + l)s'o*' ~ S'o + S'o- In turn, we can express 
5'fc+i(s'o,S'o) = S'o^ + 9 oSki 9 o, 9 o): yielding 

So(l - Sk{9*o,9o)) = 9oi^ - k9o"~')- (3-70) 

Therefore, this analysis shows that a cusped periodic wave exists iff the constants s'o and 
E satisfy the inequalities 

E > 0, S'o > 1 (3.71a) 

E>{1- kgt-^)gf >0, 9* <9^ < I (3.71b) 

^ > (1 - kgl-^)gl >0, ^<90 <9* (3.71c) 

where s'o is given by equation fl3.69p (or, equivalently, equation 03.7011 ). The solution integral 
(I3.27P then can be expressed as 

j-go „(fc-i )/2 

± / , =d9 = i (3.72) 

Jg \/{9o - 9)A{9) 

with 0 < S' < S'o- This solution describes a wave prohle s'(0 I^at has a local maximum go at 
^ = 0 modulo L and nodal cusps at ^ = ±L/2 modulo L such that S' ^ 0 while [s''] oo, 
where 

rgo „(k-i)/2 

L = 2/ / dg>0 (3.73) 

is a size parameter. 


3.4. Solitary waves. To obtain a solitary wave solution s'(0 > 0) we need 


0 > Cl > C* : E>0 (3.74) 

E = V{go) = fr(si) (3.75) 

V\go) = 0, fr"(so) < 0, (3.76) 

W(Si) > 0 (3.77) 

1/(S') < C for 0 < S'o < S' < S'!- (3.78) 

The third condition 03.761) yields 

C, = {k + l)so" - 2go (3.79) 
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and {k + l)kgQ ^ <2, which gives 


V{ 9 ) = 9 ’'^^ - / + ( 2^0 -{k + l)go)9 

(3.80) 

and 


0< 90 <9* 

(3.81) 


where g* is the point fl3.24p at which V"{g) = 0. Then the hrst part of the second condition 
fl3.75p determines 

E = gl{l-kg^,-^) (3.82) 


in terms of go. Thus, we have 

E-V{g) = ig-gofA{g) (3.83) 

and 

A{g) = 1 - ^ = 1 _ kSk{g, go) + gdgSk{g, go)- (3.84) 

9 - go 

Next, the second part of condition fl3.75p gives A{gi) = 0. This equation implicitly deter¬ 
mines gi in terms of go, as follows. We consider A{g) = 0. Since V{g) has an inflection point 
at g* < 1, we know that g > g*. For = 1, we note S'fc(l,5'o) < Sk{l, 1) = k from properties 
(| 3.39P and ( | 3.43 | ) . Hence, we have H(l)-1 = {kg^-Sk{l,go))/{l-go) > k{g^-l)/{l-go) = 
—kSk{l, go) > k‘^, which implies H(l) > 0 . But for g ^ I, clearly we have ^(^f) ~ —g^~^ < 0 
from expression fl3.84p . This implies A{g) changes sign in the interval [l,C)o) and therefore 
has a root g = gi > 1 > g*. As a consequence, the hfth condition fl3.78p is then satis- 
hed due to A{g) > 0 for go ^ g "E gi- Moreover, from expression fl3.80l) combined with 
go < g* < 1, we see that E = V{go) = goi^ — kgo~^) > 0. Hence, the hrst condition fl3.74|) 
is satished. Finally, the fourth condition fl3.77l) reduces to Sk{gi,go) > 2/{k + 1), which 
holds due to go < g* as follows. From A{gi) = 0, we get giSk{gi,go) = gi - 5'o(l - 
and so gi{Sk{gi,go) -2/{k + 1)) = gi{k - l)/{k + 1) - 5 - 0(1 - kg^-^). Then, go < g* 
gives 1 - kgl~^ > 1 - kg*^~^ = {k - l)/{k + 1 ), which yields gi{Sk{gi,go) - 2/{k 1 )) > 

{91 - 9 o)ik - l)/{k + l) > 0 , whereby we have Skigi,go) >2/{k + 1 ). 

Therefore, the preceding analysis shows that a solitary wave exists iff the constant go 
satishes the inequality (13.811) . Because go is an asymptotic minimum of the wave prohle g{^) 
as ^ ± 00 , the solution integral fl3.27|) must be expressed as 

r9l g{k-l)/2 

Jg \9 ~ 9o)\/A[g) 

with go ^ g ^ gi such that A{gi) = 0 where 5 = 51 is a (global) maximum of the wave 
prohle. The asymptotic tail has exponential decay g —)■ go, since go is a critical point of V{g) 
but V''{go) 7 ^ 0 . 


3.5. Cusped solitary waves. This type of solitary wave solution g{^) > 0 arises when 


0 > Cl > Cl: E>0 

(3.86) 

E = V{go) 

(3.87) 

V'{go) = 0 , V''{go) < 0 

(3.88) 

V{g) < E ioi Q < g < go- 

(3.89) 
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The two conditions fl3.87p - fl3.88p yield the energy expression fl3.82p and the potential 
energy expression fl3.80p . giving again the expression fl3.84p for E — V{g). But the second 
part of the third condition fl3.88p now gives 

0 < ^70 < ^ 7 * (3.90) 

where g* is the point fl3.24p at which V"{g) = 0. The fourth condition fl3.89p holds as 
follows. We note E — V{go) = 0. Then we have —V'{g) = {go — g){{k + l)Sk{g-,go) “ 2 ) = 

{2gl-^/k){go - g){{g o/g*f~^Sk{g, go) - kg^-^) g (2^ d~Vfc)(^o - g){Sk{g,go) - %o"^) < 0 
from property fl3.39p combined with the limit fl3.43p . This implies V{g) is decreasing as g 
decreases, and hence i? — > 0 for 0 < since V{g) has no critical point g < go- 

Therefore, a cusped solitary wave exists iff the inequality fl3.90p holds. The solution 
integral fl3.27p is then given by 

± / <ig = { (3.91) 

90 (5o - (/) V 9l(c/) 

with ^ < g < go, where A{g) is given by expression fl3.84p . This solution describes a wave 
profile g{^) that has an asymptotic maximum as ^ ^ ±cx), and a nodal cusp at ^ = 0 
such that S' —)■ 0 while \g'\ oo. 

The behavior of the asymptotic tail g ^ go for large ^ differs depending on whether 
Cl = Cj or Cl > Cj. These two cases respectively correspond to V"{go) = 0 and V"{go) 7 ^ 0 . 
Consequently, the general analysis using equation fl3.29p shows that go — 9 0 will exhibit 

exponential decay when V"{go) 7^ 0 , or power decay when V"{go) = 0 . 


3.6. Limiting case of periodic waves. A special limiting case of periodic waves arises 
when go = 0. This case is given by 


from which we get 
This gives 


Cl > 0 : E = 0 

(3.92) 

^(0) = V{gi) = 0 

(3.93) 

V'{gi) ^ 0 

(3.94) 

l/( 5 f) < 0 for 0 < S' < 5-1 

(3.95) 

ill exhibit nodes at which g = go = 0. 
ns fl3.92l)-fl3.93H together vield 

0 

1 

IV 

0 

(3.96) 

v{ff) = ff{ff^ - 9 i + ffi -g)- 

(3.97) 

--{ 91 - 9 ) 9 A{g), A{g) = Sk{g,gi)-1. 

(3.98) 

, we have 

9i > 1- 

(3.99) 


This implies that V'{gi) = 771 — l)>fc — 1>0 due to fc > 1, and hence the third 

condition fl3.94ll holds. The fourth condition fl3.95p now becomes A{g) > 0 for 0 < 

We first note, by direct evaluation, A(0) = Sk{0,gi) — 1 = g^~^ — 1 > 0 due to inequality 
fl3.99p . Next, we have A'{g) = S{.{g,gi) > 0 by property fl3.4ip . and hence we conclude 
A{g) > A(0) for 0 < < ^fi. This shows that condition fl3.95p holds. 
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Therefore, from this analysis, a periodic wave with go = ^ exists iff the constant gi satisfies 
the inequality fl3.99p . The solution integral fl3.27p is given by 

f9 Qk/2-l 

±/ ,/ = ( (3.100) 

Jt> y(9i - a)A{9) 

with 0 < 5 ^ < The wave profile g{^) for this solution (I3.innp will be periodic. It will have 
a local maximum of at ^ = ±L/2 modulo L, where its wavelength is given by 

rgi Qk/2-1 

L = 2 / dg>0. (3.101) 

Jo \/(9i -9)31(9) 

At .^ = 0 modulo L, the wave will have either a cusp, a corner, or a local minimum, depending 
on gi, as indicated in Table [1] which follows from inequality fl3.96p . 


3.7. Energy and momentnm. Each solution g{^) of the hrst-order ODE (I3.22p gives a 
corresponding solution of the strain ODE fl3.12p after we revert to the physical (unsealed) 
variables g{^)- From the scaling transformation fl3.16l) - fl3.18p . we get 

9(1) = l))‘''‘-0(^^l) (3.102) 

which represents a physical travelling wave in terms of the strain variable 

v = g{x-iyt) = {l{iy/cfk{k + g(^-j^^{x - ut)^ (3.103) 

where u = const, is the wave speed. Note that these waves can be bi-directional, since there 
is no constraint on the sign of u, other than z/ 7 ^ 0 . 

All physical travelling waves carry energy and momentum, which are defined by the re¬ 
spective conservation laws fl2.35p and fl2.36p where Ut = —ug and = g. The total en¬ 
ergy and total momentum are given by the spatial integrals fl2.39p and fl2.40p on a domain 
D = (—C) 0 , cx)) in the case of solitary waves, or D = (—L/2, L/2) in the case of periodic waves 
with wavelength L > 0. In terms of g{C)i these integrals take the form 


V = CM= f (3.104) 

Jn 

g = c^ig] = - c^{\R^g^-^{g'f - dl (3.105) 

For periodic waves, both integrals will be finite, since the domain is a hnite interval. 

But for solitary waves, the domain is inhnite and g{^) has a non-zero asymptotic value 
^0 > 0 as ^ ±cx), and hence both integrals will diverge. However, we can regularize the 

energy density and momentum density to remove the divergent contribution in two different 
ways. 

Firstly, we can consider a constant solution given by ^ which satisfies the strain 

ODE fl3.12l) . where go is the asymptotic value of a solitary wave solution. The respective 
momentum and energy densities fl2.35l) and fl2.36p evaluated for the corresponding constant 
strain solution v = go are given by 


Toigo] = vgh 

TAM = y~9l + i^fYo^^- 
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(3.106) 

(3.107) 










































If we now subtract these densities from the respective densities in the integrals fl3.104p and 
fl3.105p . then we obtain regularized momentum and energy integrals 


'P= [ (3.108) 

JQ 

^ ^ L ~ ~ _ ^^+ 1 ))) di ( 3 . 109 ) 

both of which will be hnite. 

Secondly, we can instead consider a linear combination of the energy and momentum 
densities such that the resulting density goes to zero as ±cxo. This is achieved simply by 
taking the density to be T^lg] — (T4[^o]/^3[^o])73[^]. From expressions (lO.lOOp and (I3.107p . 
we then have the energy-momentum integral 




k-l 


s5“‘) - di 


(3.110) 


which will be hnite. 

A direct evaluation of any of these integrals fl3.104p . fl3.105p . (13.1081) . (13.1091) . (I3.110p 
would require that we hrst substitute the explicit form of the travelling wave solution g{^) 
and then integrate with respect to ^ over the domain fl. We can avoid these obstacles by 
expressing each integral in an equivalent form (13. 6 p which involves only the ODE for the 
travelling wave. In particular, we hrst use the scaling transformation (I3.102P to express g{^) 
in terms of the dimensionless expression g{^), and we next use the dimensionless hrst-order 
ODE (I3.22P to change the integration variable from ^ to g. This formulation relies on the 
global features of g{^) that all travelling waves are symmetric functions of and that non- 
nodal periodic waves and solitary waves have a single maximum amplitude (/max at the center 
of the domain D and a minimum amplitude g'min > 0 at the domain endpoints dfl, whereas 
cusped and peaked waves as well as nodal periodic waves have a single minimum amplitude 
5'min = 0 at the center of the domain D and a maximum amplitude ^^max > 0 at the domain 
endpoints dfl. 

Thus, the momentum and energy integrals (I3.104p - (l3.105p for all periodic waves, including 
the case of cusped and peaked periodic waves, can be expressed as 

pgmax ^(fc+3)/2 

P = f dg, (3.111) 

Jgniin V ^ 

W (3.112) 

V VB - V(g) '"V 

with 


/ yillclX. 

g = ^^(fc+3)/(fc-l)|^/c|4/(fc-l)j,2^ I g{k-l)/2^ 


a = 


\l + !)• 


(3.113) 


We remark that the physical wavelength L of a periodic wave is related to the dimensionless 
wavelength L of 5 f(.^) by 


L = 


k(k+l) 

6 


LR. 


(3.114) 


Similarly, for all solitary waves, including the case of cusped solitary waves, the regularized 
momentum and energy integrals (I3.108l) - (l3.in9p along with the energy-momentum integral 
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(3.115) 


fl3.110p can be expressed as 

Jgmin '\/E — V{g) 

rgmax t _ rfl -I- . _\ 

i = 5 ^^(‘+=»/('-i)|,,/c|-/''-‘>^^fl / g(>-i)/2(g - VE-V(g)) dg, 

Jguiiu V y R y ig) / 


(3.116) 


and 


e = ^„(M)/d-i)^^/^^4ni.-i)^2jf /■'’““ j(‘-i)/ 2 ( at ‘) _ _,/( 5 ^ 

L,. V VE - v(g) ^ '"7 

(3.117) 

All of these integrals fl3.llip - fl3.112p and fl3.115p - fl3.117p are finite and represent physical 
conserved quantities for travelling waves. 


4. Explicit solutions 


The different solution integrals fl3.27p for periodic waves fl3.57p . solitary waves fl3.85p . 
cusped periodic waves fl3.72p . cusped solitary waves fl3.9ip . and nodal waves fl3.100p can be 
evaluated in an explicit analytical form in various cases which will now be systematically 
determined. 

For a given solution integral, our aim is to hnd a suitable change of variable 


g = h'^, q>0 


(4.1) 


that will bring the integral to one of the forms 



E(h) 


dh 


or 



F(h) 


dh 


(4.2) 


with < g < (y'max, where F{h) is a rational function and P{h) is a polynomial of lowest 
possible degree d > 0. If d is at most two, then the integral can be evaluated explicitly 
in terms of elementary functions; if d is three or four, then the integral can be evaluated 
explicitly in terms of elliptic functions. 

We will now classify all of the cases for which this method produces an explicit evaluation 
of each solution integral fl3.57p . fl3.85p . fl3.72p . fl3.9ip . fl3.10Up . This classihcation depends 
sensitively on the nonlinearity exponent k > 1. 


4.1. Classification of explicit solution integrals. Applying the change of variable fl4.ip 
to the general solution integral fl3.27p . we obtain 


9 


(fc-l)/2 


dg = 


qliqik+i)/2 


-1 


'50 


Jho yjE - V{hi) 


dh. 


(4.3) 


We can bring this integral to the general form fl4.2p by multiplying the numerator and 
denominator by h"^ where q and r are chosen to make be a monomial and 

(^E — V{h^)) be a polynomial of lowest possible degree. We then have E{h) = 
and P{h) = h‘^'^{E — V{h^)) in all cases when E — V{h^) has no roots with a degeneracy 
greater than one. This encompasses all of the different types of periodic waves (cf. Tabled]). 
In the case when E — V{h^) has a root h = h^ with a degeneracy of two or three, we can 
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factorize = {h — hQYB{h) to get F{h) = /iqI and P(/i) = h‘^^B{h). 

This encompasses all of the types of solitary waves (cf. Tabled]). We note that if P 7 ^ 0 then 
the lowest power term in h‘^^{E — V{h^)) is and so in this case we can assnme 2 r = 0 
or 2r = 1 (since otherwise be factored ont to get P{h) = E — V{h‘^) when 2r is 

even or P{h) = Eh — hV (h'^) when 2r is odd, which are polynomials of lower degree). 

We will now proceed to classify all possibilities for q and r snch that P{h) has degree at 
most two or at most four. If more than one possibility works for a given value of k > 1, then 
we will choose the one that yields the lowest degree for P{h). This will yield a classihcation 
of all distinct cases for which the solution integral fl4.3p can be evaluated explicitly in terms 
of elementary functions or elliptic functions. 

Note from equation fl3.24p we have 


4.1.1. Solitary waves. We begin with solitary waves fl3.85p . The most general form fl4.2p for 
the solution integral fl4.3p in this case is given by 

m = (4,5) 

h — uq 

{h - hofP{h) = - ^ 2^+2 ( 4 _g) 


with 

p=\q{k + l)Pr-l, E = hl\l-khf-^'^), Ci = -hl{2-{k + l)hf-^^). ( 4 . 7 ) 

Note we have P > 0 and Pi < 0 due to inequality (I3.8ip . Consequently, we can assume 
r = 0 or r = 1/2. Then, F{h) will be a rational function iff p is an integer, and P{h) will be 
a polynomial of degree d < 4 iff P/i^^ + is a polynomial of degree 

at most six, where we have 2r < 2r + g < 2r + 2g < 2p + 2 due to /c > 1 and g > 0. This 


determines the following classihcation: 

2p + 2 = 4, 4, 6 , 6 , 6 , 6 (4.8a) 

2r+ 2g = 2, 3, 2, 4, 3, 5 (4.8b) 

2r + g = l, 2, 1, 2, 2, 3 (4.8c) 

2r = 0, 1, 0, 0, 1, 1 (4.8d) 


where the corresponding values of d = 2 p and fc = 2 (p + l—r)/g — 1 are given by 


d = 2, 2, 4, 4, 4, 4 (4.9) 

A; = 3, 2, 5, 2, 4, | (4.10) 

The only cases that have d < 2 are when k = 2 and k = 3, where g = 1 in both cases. 

For k = 3, we have F{h) = h/{h — ho) and P{h) = {hi — h){hi + 2ho + h), where 

hi = a/I — 2h‘^ — ho with ho < h* = l/\/6. The solution integral fl4.2p in this case yields 


where 



F(h) 

Vm 


dh 


^TT — arctan 


h + ho ^ 

PW)) 


H- , arctanhf Q{h) 


Q{h) 


2^P{h)P{ho) 

P(h) + P(ho) + (h-ho)2- 
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(4.11) 


(4.12) 






















For k = 2, we have F{h) = h/{h — ho) and P{h) = {hi — h)h, where hi = 1 — 2ho with 
ho < h* = 1/3. The solution integral fl4.2p in this case yields 



m 


dh 


^TT — arctan 


h + ho — I ^ 

yp(fc) ) 


H- , arctanhf Q{h) 


(4.13) 


Hence, using h = g since g = 1, and taking into account the necessary and sufficient 
condition fl3.8ip on go, we obtain the following result. 


Proposition 1. (i) The solitary wave solution integral fl3.85p has an explicit evaluation in 
terms of elementary functions when (and only when) k = 2,3. These solutions are given by 
algebraic equations 


k = 2 : 7:71 — arctan 


9 + 9o 


+ 


9o 


arctanhl Q{g )) = ±.^, 


P{9) = {9i- 9 ) 9 , = 1-2^0, ^o</ = l/3 


(4.14) 


and 
k = 3 


1 


71 — arctan 


+ 


9o 


( 9 P 9q\ 




: arctanhl Q( 5 f)) = ±,^, 


l’(9) = ( 91 -9)(9 + 9i+ 2gi)), 9i = \/l - 29 o - 90 , 90 < 9 ’= VVfi 


(4.15) 


where Q{g) is given by expression fl4.12p . 

(a) The solitary wave solution integral fl3.85p has an explicit evaluation in terms of elliptic 
functions when (and only when) k = |,4,5. 


(i) k=2 





(ii) k=3 





Figure 7. Exact solitary wave solutions. 
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4.1.2. Cusped solitary waves. The classification for cusped solitary waves fl3.9ip is the same, 
apart from a change in the sign of the denominator of F{h) and a change in the integration 
domain. This leads to the following result. 


Proposition 2. (i) The cusped solitary wave solution integral fl3.9ip has an explicit evalu¬ 
ation in terms of elementary functions when (and only when) k = 2,3. These solutions are 
given by algebraic eguations 


k = 2 -. arctan 


{g + 9o) 


9o 


In 




•/m 

P{9) = {9i - 9 ) 9 , ^1 = 1- 2^0, 9o<9* = 1/3 


lyr -|- _ 

^ \/P{.9o) V(5'0 - 5')(1 - 25 ( 0 ) 


= ±e, 


(4.16) 


and 


k = 3 


arctan 


9o 


^/^-{9 + 9o)VW) \ 

9o{9 + S'o) + V P{9)P{^) / \/P{^ \{9o ~ 9)P{^) 


+ 


9o 


In 


9oR{9) 


= 


P{9) = { 91 -9){9 + 9i + ‘^9o), ^1 = ^/l - 2^0 - ^ 0 , go<g* = l/V6 


(4.17) 


where 


R{g) = (/^ + + ig- go?. (4.18) 

(a) The cusped solitary wave solution integral (I3.85p has an explicit evaluation in terms of 
elliptic functions when (and only when) k = |,4, 5. 


(i) k=2 





(ii) k=3 



Figure 8 . Exact cusped solitary wave solutions. 


4.1.3. Periodic waves. We next look at periodic waves fl3.57l) . In this case the most general 
form (14. 2 p for the solution integral (14. 3 p is given by 


F{h) = qhP 

P{h) = - h^P+^ 
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(4.19) 

(4.20) 
































with 


p=ig(fc + l)+r-l, E = hlhl{l-Sk{h\,hl)), Ci = Sk+i{hl,hl))-h\-hl. (4.21) 

Here we can have = 0 or Ci = 0 but not both E = Ci = 0 (since this leads to (/o = 1 
or 5 fi = 1, both of which violate inequality fl3.55p b Thus there will be three subcases to 
consider: E = 0, Ci 0] E ^ 0, Ci = 0] E 0, Ci 0. 

When E ^ 0 and Ci ^ 0, E{hi) will be a rational function iff p is an integer, and T’(h) will 
be a polynomial of degree d < 4 iff 2r < 2r + g < 2r + 2g < 2p + 2 < 4 are integers, where 
we can assume r = 0orr = l/2. This determines the following classihcation: 


2p + 2 = 4, 4 (4.22a) 

2r + 2g = 2, 3 (4.22b) 

2r + g = 1, 2 (4.22c) 

2r = 0, 1 (4.22d) 

where the corresponding values of d = 2 p + 2 and k = 2{p + 1 — r)/q — 1 are given by 

d = 4, 4 (4.23) 

fc = 3, 2 (4.24) 


There are no cases that have d < 2. 

When E = {] and Ci ^ 0, E(h) will be a rational function iff p is an integer, and P{h) will 
be a polynomial of degree d < 4 iff 2r + g < 2r + 2g < 2p + 2 < 4 are integers. Note we can 
consequently assume 2r + g = 0 or 2r + g = 1. This determines the following classihcation: 


2p + 2 = 2, 4, 4, 4, 4, 4 (4.25a) 

2r + 2g = l, 1, 2, 3, 2, 3 (4.25b) 

2r + g = 0, 0, 0, 0, 1, 1 (4.25c) 

where the corresponding values of d = 2p + 2 and fc = 2(p + l — r)/g — 1 are given by 

d = 2, 4, 4, 4, 4, 4 (4.26) 

fc = 2, 4, 2, |, 3, I (4.27) 

The only case that has d < 2 is when k = 2, where g = 1. Then we have E{h) = 1 and 


P{h) = {hi — h){h — ho) where hi = 1 — ho- The solution integral (14. 2 h in this case yields 



F{h) 


dh = arctan 



+ l^r. 


(4.28) 


When E ^ 0 and Ci = 0, E{h) will be a rational function iff p is an integer, and P{h) 
will be a polynomial of degree d < 4 iff 2r < 2r + 2g < 2p + 2 < 4 are integers. Note we can 
assume r = 0orr = l/2. This determines the following classihcation: 


2p + 2 = 2, 4, 4, 4, 4, 

2r + 2g = l, 1, 2, 3, 2, 

2 r = 0 , 0 , 0 , 0 , 1 , 
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(4.29a) 

(4.29b) 

(4.29c) 










where the corresponding values of d = 2 p + 2 and k = 2{p + 1 — r)/q — 1 are given by 

d = 2, 4, 4, 4, 4, 4 (4.30) 

fc = 3, 7, 3, f, 5, 2 (4.31) 


The only case that has d < 2 is when k = 3, where q = 1/2. Then we have F{h) = 1/2 and 
P{h) = {hi — h){h — ho) where ho = 1 — hi. In this case, the solution integral fl4.2p yields 


'ho 


^ F{h) „ 1 

— , ah = arctan 

s/m 


(izh) 

\s/mJ 


+ h- 


(4.32) 


Hence, using h = g when k = 2 since q = 1, and h = when k = 3 since g = 1/2, we 
obtain the following result after taking into account the necessary and sufficient conditions 
(I3.53p . (I3.55p . (I3.56P on gi and go. 


Proposition 3. (i) The periodic wave solution integral (I3.57p has an explicit evaluation in 
terms of elementary functions when (and only when) k = 2 with E = 0 and Ci ^ 0, and 
k = 3 with Cl = 0 and E ^ D. These solutions are given by the explicit expressions 

k = 2 : g = i(^l + (25(1 - l)cos(O), 

L = 271 


i>gi> 1/2 


(4.33) 


and 


k = 3 : 


9 = ^ 1(1 + (2c/i - l)cos(20), 

L = TT 


l>gi> l /^/2 


(4.34) 


where L denotes the wavelength. 

(a) The periodic wave solution integral fl3.57p has an explicit evaluation in terms of elliptic 
functions when (and only when) k = 2,3 with E ^ 0 and Ci 7 ^ 0; /c = |, 2, 5, 7 with E ^ 0 
and Cl = D; k = |, |, 3,4 with E = 0 and Ci 7 ^ 0. 


4.1.4. Cusped periodic waves. We continue with cusped periodic waves (I3.72p . The most 
general form fl4.2p for the solution integral fl4.3p is given by 

E{h) = qhh (4.35) 

P{h) = Eh^^ - Cih^^+^ + (4.36) 

with 

p = l{2r + q{k + 1)) - 1, Ci = CV + ^o-C, ^>0- (4.37) 

Here we can have Ci = 0, which gives two subcases to consider: Ci = 0; Ci 7 ^ 0. 

The classihcation of these two cases is the same as the corresponding two cases Ci = 0 
and Cl 7 ^ 0 for periodic waves with E ^ 0. This leads to the following result. 


Proposition 4. (i) The cusped periodic wave solution integral 03.721) has an explicit eval¬ 
uation in terms of elementary functions when (and only when) k = 3 with Ci = 0. In this 
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(i) CiT^O, E=0, k=2 (ii) Ci=0, Ey^O, k=3 



Figure 9. Exact nonlinear periodic wave solutions. 


case the solution is given by the algebraic eguation 

k = 3: g= \k(^ + (2^o “ 1) cos(2^ - 2L[i + ^/LJ)V go > 1, 

^ ^ (4.38) 

2^0 - 1 / 

where L denotes the wavelength, and [a:] denotes the floor function. 

(ii) The cusped periodic wave solution integral fl3.57p has an explicit evaluation in terms of 
elliptic functions when (and only when) k = 2,3 with Ci 0, and k = 2, 5, 7 with C\ = 0. 


L = ivT + arctan 


4.1.5. Nodal waves. Finally, we look at nodal waves (Id.lOOh . The most general form (14. 2 h 
for the solution integral 04.31) splits into two different subcases: Ci = 0 and Ci 0. 

When Cl 7 ^ 0, we have 

Fill) = qhF 
P(h) = Cih^" + 

with 

p = \qk + r — 1, Cl = hlf — h\. (4.41) 

Then F{h) will be a rational function iff p is an integer, and P{h) will be a polynomial of 


degree (i<4iff2r<2r + g<2p + 2<4 are integers. ] 
This determines the following classihcation: 

Tote 

we can assume r 

= 0 or r = 1/2. 

2p + 2 = 2, 4, 

4, 

4, 

4, 

4 

(4.42a) 

2 r + g = 1, 1, 

2 , 

3, 

2 , 

3 

(4.42b) 

2 r = 0, 0, 

0 , 

0 , 

1 , 

1 

(4.42c) 


(4.39) 

(4.40) 
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(i) Ci=0, k=3 (ii) Ci=0, k=3 



Figure 10. Exact cusped periodic wave solutions. 


where the corresponding values of d = 2p + 2 and k = 2{p + 1 — r)/q are given by 

d = 2, 4, 4, 4, 4, 4 

k = 2, 4, 2, I, 3, I 


(4.43) 

(4.44) 


The only case that has d < 2 is when k = 2, where q = 1. This yields F{h) = 1 and 
P{h) = {hi — h){h + hi — 1), from which the solution integral fl4.2p is given by 


m 

Jo l/W) 

When Cl = 0, we have 


.. 

V 2VP(/!)P(0)-(ft-i) ) 


(4.45) 


with 


F{h) = qhP 
P{h) = 1 - h^P+^ 

p= h{k-l) - 1. 


(4.46) 

(4.47) 

(4.48) 

Then F{h) will be a rational function iff p is an integer, and P{h) will be a polynomial of 
degree d<4iff0<2p + 2<4 are integers. Since p determines q = 2{p + l)/{k — 1), and 
q does appear in the powers of the terms in F{h) and P{h), we can assume p = 0 without 
loss of generality. This determines q = 2/{k — 1) and d = 2, where fc > 1 is arbitrary. Then 
we have F{h) = and P{h) = 1 — h‘^. Hence, the solution integral (14. 2 p is given by 


m 


dh = 


{k-i) 


arcsin(h). 


(4.49) 


Hence, we obtain the following result, after taking into account the necessary and sufficient 
condition fl3.99p on gi. 
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Proposition 5. (i) The periodic nodal wave solution integral fl3.100p has an explicit evalu¬ 
ation in terms of elementary functions when (and only when) k = 2 with Ci 7^ 0, and k > 1 
arbitrary with Ci = t). These solutions are given by the explicit expressions 


k = 2: g = \ - |cos(^ - L[\ + i/L\) + - 1) sin(|^ - L[\ + i/L\ |), 

1 


2 2 

L = TT + 2 arctan 


2^ S'! (s'! — 1) 


and 


k>l-. ^ = sin(i(fc-l)|^-L[i + ^/Lj|) 

27r 


2/(fc-l) 


L = 


k-l 


9i > 1, 

(4.50) 

(4.51) 


where L denotes the wavelength, and [xj denotes the floor function. 

(a) The periodic nodal wave solution integral (lO.inOh with Ci 7 ^ 0 has an explicit evaluation 
in terms of elliptic functions when (and only when) A; = |, |, 3,4. 


4.2. Properties of explicit solutions. For each of the explicit solutions derived in Propo¬ 
sitions □ toO we now exhibit their total energy and total momentum given by the integrals 
fl3.111l) - fl3.112p in the case of periodic waves and fl3.115p - fl3.117l) in the case of solitary waves. 
These integrals can be evaluated explicitly, by the same changes of variable used to evaluate 
the solution integrals. 


4.2.1. Solitary waves. From Proposition [T], we have two explicit solitary wave solutions fl4.14p 
and fl4.15p . The momentum fl3.115l) . energy fl3.116|) . and energy-momentum fl3.117p for the 
hrst solution fl4.14p are respectively given by 

r = {9/2){iy/cfliyR(^{g^ - 2go + |) + arcsin Q + | V^o(l - 3^o)^ , (4.52) 

g = {9/A){u/cfliy^R(^{l - 2go){gl - |c/o + |) + arcsin Q 

+ (fi'o ~ ffi'o + i) V 5'o(l ~ S^fo)^ , 
g = {9lfl){ylcfv^R[l - 1 ^ 0 ) (^(1 - 2 ^ 0 )^ + arcsin Q 

+ ~ tho) \/ 5'o(l “ Sf/o)^ • 


(4.54) 


(4.53) 
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(i) Ci^O, k=2 


(ii) Ci=0 



(iii) Ci=0 (iv) Ci=0 



Figure 11. Exact periodic nodal wave solutions. 

For the second solution fId.lSp . the momentum fl3.115p . energy fl3.116p . and energy- 
momentum fl3.117p are respectively given by 

V = 3\/2(z //cfvRil - 2gl) - arcsin ^ , (4.55) 

8 = (9/4)\/2(i//c)V2i?(l - 2gl) (^(l + ^^Itt - arcsin ^ - |c/o-\/l - , 

(4.56) 

£ = (3/4)\/2(i//c)V2i?(l - 2gl) - arcsin - 2gQsJl- 6g^^ . (4.57) 

The preceding expressions are plotted as functions of go in Fig. [121 
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4.2.2. Cusped solitary waves. From Proposition [2l we have two explicit cusped solitary wave 
solutions fl4.16l) and fl4.17p . The momentum fl3.115p . energy fl3.116p . and energy-momentum 
(13.11711 for the hrst solution fl4.16p are respectively given by 


-V = (9/2)(z//c)Vi?|^(5(2 - 2gQ + |) (^^ti - arcsin 


1 - 4:90 
1 ~ ‘^Qo 


V5'o(l - 2 >gQ) 


-E = (9/4)(i//c)Vi?|^(l - 2go){gl - |^o + f) - arcsin Q 
~ 2( 5^0 — fs'o + I) \/ 5'o(l — , 

-£ = (9/8)(z//c)V2i?(l - 1 ^ 0 ) (^(1 - 2 ^ 0 )^ Itt- arcsin (^^3^)) 
~ 2(1 — fs'o) V 5 'o(l ~ • 


(4.58) 

(4.59) 

(4.60) 


Note that, compared to the case of solitary waves, the overall negative sign in each of 
these expressions is due to the domain of integration being 0 < h < hmax, rather than 
hmin El h El hmax} 

Similarly, the momentum fl3.115p . energy fl3.116p . and energy-momentum fl3.117p for the 
second solution fl4.17p are respectively given by 


-V = 3V2{u/cfuR(^{l 


-goJl-3g^ , 



— arcsin 




-E = 3V2{u/cfu^R(^{gQ + g^- 



9o 

yi - 2g§ 


— arcsin 


‘^9o 

^1-2(7 



+ £/o(^i(l - t9o) \/l - 3^0 - 1(1 - 2c/o) \/ 1 - 6c/o^ ^ , 

-E = (3/4)\/2(z//c)V2i?(l - 2^2) arcsin (^-j==^ - arcsin (^-jJ== 

+ 9o (^2^1 - 69 ^ - 3^1 - 3^^^ ^ . 


(4.61) 


(4.62) 


(4.63) 


4.2.3. Periodic waves. From Proposition [3l we have two explicit periodic wave solutions 
fl4.33p and fl4.34p . The momentum fl3.11ip and energy fl3.112p are given by 

P = {9/2){u/c)\uR{gl -91 + 1 ), 

E = {9/2){iy/c)\iy^R{gl - + |) 
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(4.64) 

(4.65) 





















































for the first solution, and likewise 


(4.66) 

(4.67) 


V = (3/\/2)(z//c)^7rz/i?, 

£ = (9\/2/8) (z//c)^7rz/^i? 

for the second solution. 

The preceding expressions as functions of gi are illustrated in Fig. [T^l 


4.2.4. Cusped periodic waves. From PropositionHJ we have one cusped periodic wave solution 
fl4.38p . Its momentum (13.1111) and energy fl3.112p are given by 

V= (3/\/2)(t//c)^z/i?|^|7r + arcsin i ) + ‘^\/doido “ 1)^ > (4-68) 

£= {9V2/8){u/cfu^R(^^7i + arcsin i ) - 1)^ • (4-69) 

4.2.5. Nodal periodic waves. From Proposition [5l we have two explicit nodal periodic wave 
solutions fl4.50p and fl4.5ip . The momentum fl3.11ip and energy fl3.112p for the hrst solution 
are given by 

V = {9/2){v/cYuR(^{gl - + |) + arcsin + W9ii9i - 1)^ , (4-70) 

£ = {9/2){u/cyu^R(^{gl - + |) + arcsin (^2 g^- i ) 

+ li9i — 9i + 9 i{9i ~ 1)^ • 

Similarly, the momentum fl3.11ip and energy fl3.112p for the second solution are given by 
V = - l)-i^fc(fcTl)72^"^'^^^""'V(^)/F(f^), (4.72) 

£ = - 1)-Vfc(fc + (4.73) 

where F(x) denotes the gamma function. 


(4.71) 


5. Concluding Remarks 

For the highly nonlinear, fourth-order wave equation (II.ip . we have obtained all solitary 
wave solutions and all nonlinear periodic solutions. We are able to parameterize these 
solutions explicitly in terms of the asymptotic value of the wave amplitude in the case 
of solitary waves and the peak of the wave amplitude in the case of nonlinear periodic waves. 
All cases in which the solution expressions can be stated in an explicit analytic form using 
elementary functions are worked out. Altogether, this yields two explicit solitary waves, two 
explicit periodic waves, three explicit cusped waves, and two explicit nodal waves. We also 
give explicit expressions for the total energy and total momentum for each of these solutions. 

Our method for deriving the solutions has some novel aspects compared to the standard 
travelling wave ansatz. In particular, hrstly, we use the conservation laws admitted by the 
highly nonlinear, fourth-order wave equation (11. ip to reduce it directly to a separable hrst- 
order ODE for travelling waves. Secondly, we classify all possible types of travelling wave 
solutions by a modihed energy analysis argument. Thirdly, we use a systematic method to 
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Figure 12. (i), (ii): Energy and momentum for solitary waves, (iii), (iv): 
Energy and momentum for periodic waves. 


determine all cases of the nonlinearity exponent k > 1 for which each of the different trav¬ 
elling wave solutions can be evaluated in an explicit analytical form in terms of elementary 
functions. Fourthly, we show how to use the ODE along with its hrst integrals to obtain 
the total energy and total momentum of the travelling wave solutions without the need to 
evaluate complicated integrals involving their explicit form. This method can be applied 
more broadly to any quasilinear wave equation. 
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